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Abstract 

The repertoire of problems theoretically solvable by 
a quantum computer recently expanded to include the 
approximate evaluation of knot invariants, specifically 
the Jones polynomial. The experimental implementa- 
tion of this evaluation, however, involves many known 
experimental challenges. Here we present experimental 
results for a small-scale approximate evaluation of the 
Jones Polynomial by nuclear-magnetic resonance (NMR) , 
in addition we show how to escape from the limitations 
of NMR approaches that employ pseudo pure states. 

Specifically, we use two spin 1/2 nuclei of natural abun- 
dance chloroform and apply a sequence of unitary trans- 
forms representing the Trefoil Knot, the Figure Eight 
Knot and the Borromean Rings. After measuring the 
state of the molecule in each case, we are able to esti- 
mate the value of the Jones Polynomial for each of the 
knots. 



INTRODUCTION 

The Jones polynomial [T , a great discovery in knot the- 
ory, has recently become an interesting topic for quantum 
computing. In particular, the use of quantum comput- 
ing has been discussed for approximately evaluating the 
Jones polynomial V{z) at selected values of z. For a knot 
displayed as a braid of n strands (specified in terms of 
a sequence of crossings), these are the values z of the 
form z = exp(27ri/fc) where k is an integer in the algo- 
rithm of Aharonov, Jones and Landau (AJL) [5]. In [3111] 
a quantum algorithm is given by Kauffman and Lomo- 
nanco (KL) for three-strand braids that can be used to 
evaluate the Jones polynomial at a continous range of an- 
gles. Most of the computational cost of the approximate 
evaluation is the estimation of the trace of a unitary ma- 
trix. The method of estimation described for the AJL 
algorithm and the KL algorithm requires that the quan- 
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tum computer separately obtain an estimate of each of 
the diagonal elements of the unitary matrix; then these 
estimates are summed to yield an estimate of the trace. 

The next section reviews the relation of the Jones poly- 
nomial for a braid to a unitary transformation composed 
of factors that correspond to braid crossings, so that the 
problem of evaluating the Jones polynomial reduces to 
the problem of evaluating traces of unitary matrix. 

An exposition of how the KL algorithm (which we use 
in this paper) can be regareded as a special case of a 
generalization of the AJL algorithm is presented after 
that. In this sense this paper and its sequels will be 
about experimental implementation of both the KL and 
the AJL quantum algorithms for computing the Jones 
polynomial. 

Following this, we present the method whereby an ide- 
alized NMR quantum computer [5l |6] can evaluate the 
trace of unitary matrix written as a product of factors 
all at once, that is, without having to evaluate diagonal 
elements of the unitary matrix separately. 

Experimental results for the evaluation of cases of the 
2-by-2 matrix, and hence of the Jones polynomial for a 
braid of 3 strands, by use of nuclear magnetic resonance 
(NMR) is the subject of the last section. 



THE JONES POLYNOMIAL AND UNITARY 
MATRICES 

The Jones polynomial [T] was a great discovery in knot 
theory. It marked the beginning of a significant rela- 
tionship between knot theory and statisical mechanics, 
particularly through the relationship of the polynomial 
with the Temperley-Lieb algebra, and through the ex- 
plicit bracket state sum model 13 [SI IB [IHl [II]- From 
the topological side the Jones polynonmial is striking be- 
cause it can detect the difference between many knots 
and their mirror images. 

The general algorithm to find the Jones polynomial is 
in the jJP complexity class, and so this is an algorithm 
worth understanding in the context of quantum compu- 
tation. 
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The key idea behind the present quantum algorithms 
to compute the Jones polynomial is to use unitary rep- 
resentations of the braid group derived from Temperley- 
Lieb algebra representations that take the form 

p{ai) = AI + A-^U^ 

where Ui is a standard generator of the Artin braid group, 
A is a complex number of unit length, and Ui is a sym- 
metric real matrix that is part of a representation of 
the Temperley-Lieb algebra. For more details about this 
strategy and the background information about the Jones 
polynomial, the bracket model for the Jones polynomial 
and the Temperley-Lieb algebra the reader may wish to 
consult [llEllSlllliaiHlElIlOliniinilTaj. In the follow- 
ing mathematical description, we have given a minimal 
exposition of the structure of such representations. 



Two Projectors and a Unitary Representation of the 
Three Strand Braid Group 

It is useful to think of the Tcmperlcy Lieb algebra as 
generated by projections Ci = Ui/ 6 so that = and 
e^Ciiiei = TCi where r = 
K-J|>1- 



and ei and Cj commute for 



With this in mind, consider elementary projectors 
e = \A)(A\ and / = |-B>(-B|. We assume that (A\A) = 
{B\B) = 1 so that = e and p = f. Now note that 



efe = \A){A\B){B\A){A\ = {A\B){B\A)e = re 



Thus 



efe — re 



where t = {A\B){B\A). 

This algebra of two projectors is the simplest instance 
of a representation of the Temperley Lieb algebra. In 
particular, this means that a representation of the three- 
strand braid group is naturally associated with the alge- 
bra of two projectors. 

Quite specifically if we let {A\ = {a,b) and \A) — 
{a,b)^ the transpose of this row vector, then 



Note also that eie2 — 
We define 



ah 




and 6261 — 



a' 
ab 



Ui = 5ei 

for i = 1, 2 with a? = 6^"^ . Then we have , for i 
Uf = SU, , U1U2U1 = C/i, U2U1U2 = U2. 
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Thus we have a representation of the Temperley-Lieb al- 
gebra on three strands. See [lU] for a discussion of the 
properties of the Temperley-Lieb algebra. 
Note also that we have 



while 



trace{Ui) — trace{U2) = S, 



trace{UiU2) = trace{U2Ui) = 1 



where trace denotes the usual matrix trace. We will use 
these results on the traces of these matrices in Section . 

Now we return to the matrix parameters: Since a^ -\- 
fe2 — 1 tJijg means that 5~'^ + b'^ — 1 whence 6^ = 1 — 5^^. 
Therefore b is real when 5"^ is greater than or equal to 1 . 

We are interested in the case where 5 — —A^ — A^"^ and 
A is a unit complex number. Under these circumstances 
the braid group representation 

p{ai) ^ AI + A-^Ur 

will be unitary whenever Ui is a real symmetric matrix. 
Thus we will obtain a unitary representation of the three- 
strand braid group B^ when 5^ >\. 
For any A with d = —A? — A^^ these formulas define a 
representation of the braid group. With A — exp{i9), we 
have d = —2cos{29). We find a specific range of angles 9 
in the following disjoint union of angular intervals 

e e [0,7r/6]U[7r/3,27r/3]U[57r/6,77r/6]U[47r/3,57r/3]U[ll7r/6,27r] 

that give unitary representations of the three-strand braid 
group. Thus a specialization of a more general represen- 
tion of the braid group gives rise to a continuous family 
of unitary representations of the braid group. 



ab 
ab &2 



is a standard projector matrix when a^ 
obtain a specific representation. 



1. To 



let ei — 



1 




and 62 — 



a^ ab 
ab 62 



It is easy to check that 616261 = a^ei and that 626162 



0^62. 



A Quantum Algorithm for the Jones Polynomial on 
Three Strand Braids 



We gave above an example of a unitary representation 
of the three-strand braid group. In fact, we can use this 
representation to compute the Jones polynomial for clo- 
sures of 3-braids, and therefore this representation pro- 
vides a test case for the corresponding quantum compu- 
tation. We now analyse this case by first making explicit 
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how the bracket polynomial is computed from this rep- 
resentation. This unitary representation and its applica- 
tion to a quantum algorithm first appeard in [3]. When 
coupled with the Hadamard test, this algorithm gets val- 
ues for the Jones polynomial in polynomial time in the 
same way as the AJL algorithm [5]. It remains to be seen 
how fast these algorithms are in principle when asked to 
compute the polynomial itself rather than certain spe- 
cializations of it. 

First recall that the representation depends on two ma- 
trices Ui and U2 with 
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and U2 = 



The representation is given on the two braid generators 

by 



and 



p{ai) = AI + A-^Ui 



p{a2) = AI + A-^U2 



(1) 



(2) 



for any A with 5 = —A^ — A 2, and with A = exp{i6), 
then i5 = —2cos{29). We get the specific range of angles 
9 € [0,7r/6] U [7r/3,27r/3] U [57r/6, 77r/6] U [47r/3, Stt/S] U 
[IItt/G, 27r] that give unitary representations of the three- 
strand braid group. 

Note that tr{Ui) = tr{U2) = 5 while tr{UiU2) = 
tr{U2Ui) = 1. If b is any braid, let 7(6) denote the sum 
of the exponents in the braid word that expresses h. For 
h a three-strand braid, it follows that 

where / is the 2x2 identity matrix and r(6) is a sum of 
products in the Temperlcy Lieb algebra involving Ui and 
U2- Since the Temper ley Lieb algebra in this dimension 
is generated by /,C/i, U2, U1U2 and U2U1, it follows that 

(b) = A^('')<52 + ir(T(6)) 

where b denotes the standard braid closure of b, and 
the sharp brackets denote the bracket polynomial as de- 
scribed in previous sections. From this we see at once 
that 

(b)=trip{b))+A'^''\5^-2). 

It follows from this calculation that the question of com- 
puting the bracket polynomial for the closure of the 
three-strand braid b is mathematically equivalent to the 
problem of computing the trace of the matrix p{b). 

The matrix in question is a product of unitary matri- 
ces, the quantum gates that we have associated with the 
braids ci and a2- The entries of the matrix p{b) are the 
results of preparation and detection for the two dimen- 
sional basis of qubits for our machine: 



Given that the computer is prepared in \j), the proba- 
bility of observing it in state \i) is equal to 
Thus we can, by running the quantum computation re- 
peatedly, estimate the absolute squares of the entries of 
the matrix p(b) . This will not yield the complex phase in- 
formation that is needed for either the trace of the matrix 
or the absolute value of that trace. 

However, we do know how to write a quantum al- 
gorithm to compute the trace of a unitary matrix (via 
the Hadamard test). Since p{b) is unitary, we can use 
this approach to approximate the trace of p{b). This 
yields a quantum algorithim for the Jones polynomial 
for three-stand braids (evaluated at points A such that 
the representation is unitary). Knowing tr(p(b)) from 
the quantum computation, we then have the formula for 
the bracket, as above, 

(b) ^ trace{p{b)) + A'^^\5^ - 2). 

Then the normalized polynomial, invariant under all 
three Reidemeister moves is given by 

/(6) = (-A3)-^W(S). 

Finally the Jones polynomial in its usual form is given 
by the formula 

v{b){t) = fm-^''). 

Thus we conclude that our quantum computer can ap- 
proximate values of the Jones polynomial. 



ON THE RELATIONSHIP WITH THE AJL 
ALGORITHM 

Here is how the KL (Kauffman-Lomonaco) algorithm 
described in the previous section becomes a special case 
of a generalization of the AJL algorithm: Here we use 
notation from the AJL paper. In that paper, the genera- 
tors Ui (in our previous notation) for the Temperley-Lieb 
algebra, are denoted by Ei. 

Let Lk — Afc ~ sin{kO). For the time being is an 
arbitrary angle. Let A = iexp{i6/2) so that d = —A^ — 
= 2cos{0). 

We need to choose 9 so that sin{k9) is non-negative 
for the range of k's we use (these depend on the choice of 
line graph as in AJL). And we insist that sin{k9) is non- 
zero except for fc = 0. Then it follows from trigonometry 
that {Lk-i + Lk+i)/Lk — d for all k. 

Recall that the representation of the Temperley-Lieb 
algebra in AJL is given in terms of Ei such that Ef = dEi 
and the Ei satisfy the Temperley-Lieb relations. Each 
Ei acts non-trivially at the i and i + I places in the 
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bit-string basis for the space and each Ei is based upon 
La-i, La, La+i where a = z{i) is the endpoint of a walk 
described by the bitstring using only first {i — 1) bits. 
Bitstrings represent walks on a line graph. Thus 1011 
represents the walk Right, Left, Right, Right ending at 
node number 3 in 

1 2 3 4. 

For p = 1011, z(l) = l,z{2) = 2,z(3) = l,z(4) = 
l,z(5) = 3. 
More precisely, if we let 

\via)) = [^La-l/La, ^/La+l/Laf 

(i.e. this is a column vector. T denotes transpose.) Then 

E., = \v{zm(v{zm 

Here it is understood that this refers to the action on the 
bitstrings 

01 

and 

10 

obtained from the given bitstring by modifying the i and 
i+1 places. The basis order is 01 before 10. Conceptually, 
this is a useful description, but it also helps to have the 
specific formulas laid out. 

Now look at the special case of a line graph with three 
nodes and two edges: 

1 2 3. 

The only admissible binary sequences are |110) and |101), 
so the space corresponding to this graph is two dimen- 
sional, and it is acted on by Ei with z{l) = 1 in both 
cases (the empty walk terminates in the first node) and 
E2 with z(2) = 2 for 1 110) and z(2) = 2 for |101). Then 
we have 

-BilllO) = 0,£;i|101) (i|101), 

E2\xyz) = \v)(v\xyz) 
(xyz = 101 or 110) where v = {^/ijd, ^d-l/d)^. 

If one compares this two dimensional representation 
of the three strand Temperley - Lieb algebra and the 
corresponding braid group representation, with the rep- 
resentation Kauffman and Lomonaco use in their paper, 
it is clear that it is the same (up to the convenient re- 
placement of ^ = exp{i9) hy A = iexp{id/2)). The trace 
formula of AJL is a variation of the trace formula that 
Kauffman and Lomonaco use. Note that the AJL algo- 
rithm as formulated in [5] does not use the continuous 
range of angles that are available to the KL algorithm. 
In the sequel to this paper and in a separate paper on 
the mathematics, we shall show how the entire AJL al- 
gorithm generalizes to continuous angular ranges. 



THEORY OF AN NMR SPECTROMETER USED 
AS A QUANTUM COMPUTER 

By convention, a quantum computer as conceived in 
theory is assumed to yield an outcome associated with 
a quantum measurement of some (possibly mixed) quan- 
tum state. In contrast, NMR machines implement a re- 
stricted version of an Expectation- Value Quantum Com- 
puter (EVQC), which in place of an outcome yields, to 
some finite precision, the expectation value for a mea- 
surement of a (again, possibly mixed) quantum state [5]. 
Reflecting facts of NMR spectrometers, an NMR Quan- 
tum Computer (NMRQC) implements only the special 
measurement operators discussed in [14] , and these mea- 
surement operators all have zero trace. 

Here are the details. For a Hermitian measurement 
operator M applied to a density matrix p, the EVQC of 
precision e yields a value x such that 

|a;-Tr(A/p)| < eA(M), (3) 

where A(M) is the difference between the minimum and 
the maximum eigenvalue of the measurement operator 
M, which is just the possible range to the trace as p varies 
over all possible density matrices. (The factor A(M) 
makes limitations of resolution immune to the mere an- 
alytic trick of multiplying the measurement operator by 
a constant.) 

The measurement operators of main interest for the 
algorithm by which we estimate the trace of a unitary 
operator are Ii^ and Iiy, shortly to be defined. 

Thermal Equilibrium and initial state preparation 

To first order, the initial thermal state density operator 
of an ensemble (very large number) of quantum systems 
with n + I qubits each |15j is given by 

^ n+l 

m« ^(l-^tti/z.) (4) 

1=1 

with ai = ^, 

= ^1 «) ■ ■ ■ ® 1 «) cTz «) 1 «) ■ • ■ «) 1, 

(where the Pauli matrix ct^ appears as the l^^ term in the 
product), uji is the resonance frequency of qubit Z, k is 
Boltzmann's constant, T is temperature and N = 2"+^. 

The initial density operator required for our algorithm 
is given by 

Po = ^(1-ai/iz) (5) 
which can be prepared from pth by a variety of methods 
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Algorithm to estimate the trace of U 

The method presented here is based on the algorithm 
that first appeared in [I7|. As mentioned above we as- 
sume that U is given in the form of local operations on 
n qubits. Given a program for U, Barenco et al. |18j 
describe a procedure to construct a program or local op- 
erations for the operator controUed-U, cU. cU operates 
on n -I- 1 qubits, does not affect the first qubit, applies 
U on the remaining n qubits if the first qubit is |1) and 
does nothing otherwise: 

cu\im = (i0c/)|i)iv) = ii)t/i^) 
cu\om = io)i^). 

In block matrix form, cU is given by: 

1 



U 



We now describe our algorithm: 
Step 1: Prepare the density operator: 

1 , . 1 ai / 1 

^ ' N 2N \1 

where 



(6) 



hx = 2^a; (8) 1 and cr^ 



Step 2: Apply cU to pi: 



1 

1 



p,^cUp,cU^^-l-—l Q 



Step 3: Measure {Iix + Hiy) to estimate 

trace{{hx + ihy) P2) = —trace{U), 



(7) 



(8) 



where 



Iiy — -ay 1 and ay = 



-i 

1 



By Eq. (|3| the result of this measurement is a complex 
number, z, such that 



|Re(z) 



N 



Re[trace{U)]\ < e 



|Im(z) - -j^lin[trace{U)]\ < e. 



It follows that the measurement result satisfies 
\z - trace{U)\ < V2Ne/ai. 



(9) 



(10) 



EXAMPLE KNOTS AND EXPERIMENTAL 
RESULTS 

Experimental results for 3 knots on three strands were 
obtained using the methods outlined above. Specifically, 
we present results for the Trefoil Knot, the Figure Eight 
Knot, the Borromean Rings. A 2-spin system (details 
of the molecule and pulse sequences are given later) was 
used, the initial state given by the density operator pro- 
portional to was prepared and a reference spectrum 
was then collected. This was followed by application of 
a controUed-unitary operator corresponding to and rep- 
resenting each knot separately found from the represen- 
tation 

si=p{ai) and S2 = p{a2). 

Measurement of the expectation value of Iix + iliy after 
applying the controUed-unitary operator yields the trace 
of the unitary operator representing the knot and thus 
the estimate of the Jones Polynomial for each knot. 

For each of the three knots, the Jones Polynomial was 
estimated at the complex numbers e*^ for all 9 in the 
range < < tt/6 at single degree increments (31 val- 
ues). Comparison to the theoretical values shows ex- 
cellent correspondence with experimental observations. 
Furthermore, the Jones Polynomial itself for each of these 
knots can be constructed from the experimental results. 



Experimental setup and molecule 

All experiments were performed on a Bruker Avance 
DMX 750 NMR spectrometer, equipped with a TXI 5mm 
probe head with XYZ gradients. The sample was a 9:1 
mixture of chloroform and deuterated acetone. It natu- 
rally contained about 1% of ^'^C-^H chloroform which 
was the active compound that represented the hard- 
ware of our NMR quantum computer. The spin systems 
Larmor frequencies were 188.6349005 MHz for ^'^C and 
750.1354275 MHz for ^H. The corresponding chemical 
shifts are 77.2 ppm and 7.235 ppm, respectively. The 
two spin-^/2 nuclei of ^^C-^H chloroform interact through 
scalar coupling. The corresponding coupling constant is 
J=209.5 Hz. The longitudinal relaxation times (Ti) and 
transversal relaxation times (T2) of both spin-^/2 nuclei 
are: Ti: 21.8 sec, T2: 0.19 sec, ^H Ti: 6.1 

sec, ^H T2: 0.48 sec. In order to suppress the signal of 
99% ^^C-^H chloroform and to prepare the initial opera- 
tor Ix, where / corresponds to ^'^C (^H will be referred 
to as S), the following preparation sequence was used in 
all experiments: the ^H spins were saturated by cw ir- 
radiation. Subsequently they were dephased by applying 
a 9.9 /Lts 90° (^H) pulse followed by a Bp gradient. This 
sequence of 90° pulse and gradient was repeated twice 
with orthogonal gradients. Subsequently the ^'^C spin 
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was excited using a 19.45 /iS 90° (^^C) pulse. This prepa- 
ration sequence was followed by the pulse sequence of 
the individual experiments (see Figure 7). Finally the 
"'^'^C signal was detected by measuring 512 points during 
452 ms. In order to improve the sensitivity, we decoupled 
all spins during the detection period by applying the 
DIPSI-2 [IH] decoupling sequence. 



CONCLUSION 

In this paper, we showed how the KL algorithm is a 
special case of a generalized AJL algorithm. Using the 
KL algorithm, we obtained a unitary representation of 
the three-strand braid group and discussed a method for 
computing the Jones polynomial using this representa- 
tion over a range of complex numbers. Next, the theory 
of an idealized NMR quantum computer was presented 
and we showed how the trace of a unitary matrix can be 
experimentally determined. Experimental realization for 
three different knots where performed where the exper- 
imental data agreed with theoretical calculations. Fu- 
ture work includes generalizing the AJL algorithm for 
any number of strands, as was done for the three-strand 
braid group in this work, and their experimental imple- 
mentations. 
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FIG. 2: Experimental results for the Trefoil Knot FIG. 4: Experimental results for the Figure Eight Knot 




FIG. 3: The Figure Eight Knot generated by the sequence FIG. 5: The Borromean Rings generated by the sequence 

-1 -1 -1 -1 -1 
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FIG. 6: Experimental results for the Borromean Rings 



FIG. 7: Pulse sequence to implement a controlled-si,2 op- 
eration (left) and its inverse (right). For si set 7 to 0. 
(To obtain a propagator of 1 © si,2 respectively its inverse, 
we have to apply a global phase factor of e^"''''''^. The 
propagator of a pulse on the second spin- 1/2 is defined as 
Uprise ■- e-*'^% with S„ ■- (l®cr^/2),e G {a,/3,7} 
and V £ {y, z} where a = O.Stt — 26, /3 = O.Svr -I- and 
7 = tan*^[cos46l/V4cos2 29- 1] + n/2 and < 6» < 7r/6. 
For the free evolution UevoiuUon ~ g-^^tJ^i^s^ ^j^j^ 21^3;, := 
2{az/2 ® crz/2). The pulse sequence is applied to the initial 
density operator Ii^ 



